Vibrational states and disorder in continuously compressed model glasses 
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We present in this paper a numerical study of the vibrational eigenvectors of a two-dimensional 
amorphous material, previously dee ply studied from the point of view of mechanical properties 
and vibrational eigen- frequencies [7Hl0l|. Attention is paid here to the connection between the 
mechanical properties of this material in term of elastic heterogeneities (EH), and how these inherent 
heterogeneous structures affect the vibrational eigenvectors and their plane waves decomposition. 
The systems are analysed for different hydrostatic pressures, and using results from previous studies, 
a deeper understanding of the boson peak scenario is obtained. The vibrational spectrum of a 
continuously densified silica glass is also studied, from which it appears that the pulsation associated 
with the boson peak follows the same pressure dependence trend than that of transverse waves with 
pulsation associated with the EH characteristic size. 

PACS numbers: 61.43.Fs; 63.50.Lm; 62.30.+d. 



I. INTRODUCTION 

Glasses develop a resistance to a macroscopic defor- 
mation, like solids, but depict a lack of structm'al order, 
like liquids. In solids, the response to such a deforma- 
tion is theoretically approached by the use of the con- 
tinuum elasticity, which involve phonons to model the 
vibrational excitations. Due to the loss of translational 
invariancc in glasses in comparison to the corresponding 
crystal, propagating phonons are difficult to describe, but 
have already been observed down to the atomic level [l| . 
It's also admitted that a common characteristic of dis- 
ordered solids, when compared to crystalline ones, lies 
in the existence of an excess of low-frequency vibrational 
states, the well known boson peak. This excess produces 
the also well-known specific heat anomaly of glasses at 
temperatures T ~ 10 K, and seems to be a signature of 
the disorder in glasses beyond the nanometer scale. 

Many attempts have been proposed to theoreti cally 
interpret the boson peak anomaly [l-lijp, HI HI, [sg , 
and among them, Shintani and Tanaka have recently 
proposed an interpretation in term of involved atomic ar- 
rangements with transverse vibrations that tend to be 
particularly strong, in two-dimensional systems, around 
favoured structures of five-fold symmetry. The frame- 
work in which the boson peak is more related with 
transverse vibrations than longitudinal ones has been 
already approached in a series of papers dealing with 
the vibrational and mechanical properties of two- and 
three-dimensional soft (Lennard- Jones) and strong (sil- 
ica) glasses [vl-fioj. In those systems, transversal motions 
of the particles give rise to large vortices of characteristic 
size inaff, when they respond to a macroscopic deforma- 
tion. This noisy displacement field has been called non- 
affine displacement field, and involve correlated displace- 
ment of atoms from S,naff ^ 20 inter-atomic distances 
in three dimensions, to £,naff 30 in two dimensions. 
In both soft and strong glasses, the application of the 
classical continuum elasticity theory is subject to strong 
limitations below such length scales Cna//- This effect is 
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particularly important when the wavelength of the vibra- 
tional excitation is less than this characteristic size ^„a/ / ■ 
The origin of the departure from classical behaviour is 
very likely related to the disorder in inter-atomic interac- 
tions (local stresses, inhomogeneities in elastic constants, 
local anisotropy) . This leads to an important decrease of 
the average shear modulus in amorphous systems com- 
pared with the corresponding crystals. Note that the 
typical size ^na// depends on the pressure: it decreases 
and saturates at high pressures for both soft and strong 
glasses. Finally, it has been also previously noticed that 
the estimate of the frequency associated with S^nojj is in 
good agreement with the boson peak position [lO[, and 
thus should encourage to consider the boson peak as a 
length marking a crossover between (i) a. regime where 
vibrations with wavelengths larger than £,naff niay be 
well described by a classical continuum theory, and (ii) 
a small wavelength regime where vibrations are strongly 
affected by the elastic heterogeneities (EH). 

In the following, such a picture will be confirmed in the 
case of a two-dimensional glass. After a summary of the 
theoretical concepts related with the continuum elastic- 
ity and an overview of the simulated system, the pressure 
dependence of the vibrational density of states (vDOS) 
and the boson peak will be approached for both strong 
and soft glasses. Wc then expect significant change of 
the vDOS upon compression, namely a decrease of boson 
peak height and its shift toward higher frequencies [l^ . 
While this kind of pressure dependence study is merely 
used, in the literature, to im pro ve the theoretical inter- 
pretations of the boson peak [l3l - [l6l [2ll . Is^ , most of the 
present work will be devoted in the effect of increasing 
pressure, on the interplay between clastic heterogeneities 
(EH) and clastic wave propagation in a two-dimensional 
glass. 



II. TECHNICAL DETAILS 

We first briefly review the considered systems studied 
in this work, much more details about them being given 
in previous works Then, the theoretical frame- 

work used in background of this work is summarized. 
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A. Simulated systems 

The two-dimensional prototype materials we consider 
are bulk-like systems contained in a square of side L with 
periodic boundary conditions. Computationally, they are 
formed by quenching very quickly a slightly polydisperse 
liquid of iV = 10000 spherical particles interacting via 
simple Lennard- Jones (LJ) pair potential, into the near- 
est energy minimum, following a fixed protocol 0, Q us- 
ing standard molecular dynamics, steepest descent and 
conjugate gradient. The resulting structures are amor- 
phous, i.e. they exhibit not a crystal-like, but a liquid- 
like order. The density of the bulk-like 2D glassy mate- 
rials in the final state arc adjusted in order to cover a 
wide range of hydrostatic pressure from P — Ocr^/e to 
P ~ 130tT^/e (in LJ units). This kind of glass mimics a 
soft glass according to Angell's classification (T^]. On the 
other hand, an amor pho us silica glass is also simulated in 
the BKS framework [l8| , and using a fixed quench proto- 
col that do not alter the final glassy state. For this 
kind of glass, systems containing between N = 24048 
and N = 64128 particles were used, and for hydrostatic 
pressures from P = GPa to P = 20 GPa. 



B. Vibrational modes from the continuous 
elasticity theory 

In the framework of the linear classical theory of elas- 
ticity, which assume that elastic bodies are homogeneous 
and isotropic, the Landau free energy can be written on 
the form: 



= Y. dk{2TT)"'S{^ - ljM) (3) 

Using the Debye's model [l^ for set of phonon propa- 
gating in a continuum medium, this vibrational density 
of states can be explicitly evaluated, which leads to the 
famous formulae: 

(4) 

In particle based simulations, it's in principle possi- 
ble to access the low temperature vibrational spectrum 
of glasses using the harmonic approximation Uxot = 
Y^iLi^Ljirij) « Y.i,j Mia,ji3Ua{ri)ui3{rj). Heucc, the 

equation of motion for an atom i, rrii ^g"2° (r J — — , 
becomes equivalent for the whole system to an Euler- 
Lagrange eigenvalue problem (mw^ — M ) .1 = on the 
plane wave basis. Therefore, the vibration frequencies 
around an equilibrium position can be exactly calculated 
by diagonalizing the dynamical matrix M, which is ex- 
pressible in terms of the first and second derivatives of 
the inter-particle interaction potentials, while the corre- 
sponding displacement fields are given by the eigenvec- 
tors. 



III. RESULTS FOR THE VIBRATIONAL 
SPECTRUM 
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where (A, fj,) are the Lame coefhcients, d the space di- 
mension, and the strain eap = {dpUa + daUp) /2 for 
the continuum displacement field u. The stress is re- 
lated to the strain using the classical relation ctq,^ = 
'^ai3 + '^TreJc^ -|- 2^eap , where a^p is the reference state 
stress. Using this formalism, the equation of motion for 

the elastic field pduUa = Yl'p=i^P'^oip{^) + be 
written in a wave formulation: 



The Eq. (0) can be solved using the Stokes decomposi- 
tion u = V4> -I- V X ■)/', with (p = Cfy'^(l) for longitudinal 
waves of speed C£ = (A + 2y) /p, and i/j = C^\7^ip for 
transverse waves of speed = ii/p. These two wave 
equations can be solved using the plane waves basis with 
wave length Ap = 2tt/ \\ k \\. For finite size systems as the 
ones we're simulating, d quantum numbers appear and 
for instance in two-dimensions, (fc^;, ky) — {2tt/L) {n,m). 
Then, the vibrational density of states (vDOS) is ex- 
pressed as: 



Using the formalism developed in the previous Sec lIIBl 
a first look on vibrational eigen-frequencies is presented 
in the following, by also varying the applied hydrostatic 
pressure of simulated systems. The three-dimensional 
glass is approached using a silica glass, and the two- 
dimensional one using a Lennard- Jones glass. 



A. Pressure dependence of the vibrational density 
of states 

On FigUfa), the vDOS from Eq.® for a silica glass 
under different hydrostatic pressure are plotted. Results 
are in agreement with numerical and experimental lit- 
erature [23, m m iH, Hi]. In the inset of this figure, 
the ratio to the Debye's prediction Eq.Q which depicts 
the boson peak anomaly is shown. It appears that the 
boson peak intensity decreases with increasing pressure 
and tends to be unmeasurable at very high pressure. The 
corresponding bos on p eak frequency i^bp also shifts to 
higher frequencies [IlSIlil. 

On the Figl^Ja), a scaling analysis for the silica glass 
is plotted using reduced units v/vbp, which may lead 
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to a collapse of the vDOS on a single master curve 
at least depending the strength of the hydrostatic pres- 
sure. We then observe for pressures less than 7 GPa, 
that the rescalcd vDOS collapse well (as also observed in 
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Fig 1. Vibrational density of states (vDOS) for (a) a sil- 
ica glass at various pressures, and (b) for a two-dimensional 
Lennard- Jones glass (LJ) at a density corresponding to a null 
hydrostatic pressure state. In the inset of both figures, the 
ratio to the Debye's prediction Eq.Q is plotted, which shows 
the excess of vibrational density of state. On the (b) part, 
the estimated contribution of longitudinal (L) and transverse 
(T) pulsations uit/l{P) = '^''^Ct / l{P) / £,naf f{P) for the char- 
acteristic size $,naff{P) of the elastic heterogeneous do- 
mains are also drawn. 



Fig 2. Vibrational density of states for the silica glass and for 
various pressures, plotted in v/vbp[P) rescaled units, where 
vbp{P) is the boson peak frequency. Both scaling lead to 
a curve collapse for pressures less than 7 GPa, above which 
the silica glass involves strong structural modifications, (a): 
vDOS are rescaled according to [2l|]. Inset: relative intensity 
of boson peak and Debye level, normalized to unity by P = 
GPa values for each quantity. The dashed line marks the 
plastic transformation threshold upon compression [SJ. (b): 
vDOS are rescaled to conserve modes [22|. Inset: rescaled 
boson peak and Debye frequency (see text) versus pressure. 



21[), but for higher pressures, this collapse is lost. This 
is related to a great change in the structure of the con- 
sidered strong glass former at such pressures [1^ [l^l ■ In 
the inset of the Figl^Ja) , the variation of the boson peak 
intensity Ibp{P) is compared to the one of the Debye 
level, while the onset of irreversible compression around 
8.6 GPa is also drawn [l^. In this inset, quantities are 
normalized to unity by their values at P = GPa. It has 
been observed experimentally [32| . that the boson peak 
shift follows the variation of the Debye level under perma- 
nently densified silicate glasses, and so should be related 
to a transformation of the elastic medium. A contrary in 
[2ll |. the relative variation between both quantities has 
been shown to increase, i.e. the boson peak shift pres- 
sure dependence is stronger than for the Debye level. The 
authors in 2l| thus conclude that the boson peak phe- 
nomena may not be only explained by a transformation 
of the elastic continuum upon compression. 



In the inset of the Figl^^a), the relative ratio between 
Ibp{P) and the Debye level decreases with the same 
trend upon compression, which should also corroborate 
the full continuum medium transformation (CMT) pic- 
ture. In fact, the situation is more complex, as shown 
on the FigEl where the CMT statement g^^" {lo)/uj'^ cx 



{M)~^ [C^^ + 2C^^) is plotted versus the 



{d{v))-' 

compression level. As already pointed out in [12|, it's 
observed a strong non-linear pressure dependence of the 
CMT prescription in this silica glass. This suggests that 
the CMT picture cannot fully, or at least solely, explain 
the boson peak nature and behaviour under compres- 
sion. A confirmation is then obtained on the FiglJl^b), 
where the vDOS arc plotted using another set of rescal- 
ing units g{v)-i'%p/v^ that takes into account the con- 
servation of modes, while keeping the adimensional pa- 
rameter V jvBP on abscissa. Again, reasonable collapse of 
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Fig 3. Variation of the continuous medium transformation 
(CMT) statement with hydrostatic pressure for a silica glass. 
A complex trend is observed, with a maximum close to the 
one measured in This points out that the CMT cannot 
solely explain the boson peak dependence with pressure. The 
dashed line is a guide for the eyes. Inset: pressure variation 
of the ratio between identical rotations of inter-tetrahedral 
units and opposite ones, after an elastic elongation has been 
applied to a silica glass. At high pressure, identical rotations 
dominate, leading to a decrease of inhomogeneous particles re- 
arrangements. Connected with the decrease of elastic hetero- 
geneities characteristic sizes, this could explain the weakness 
of the boson beak intensity, while its shift could be approached 
by the one of speed of sounds for transverse waves. 



vDOS is achieved for pressures less than 8.6 GPa. More 
intriguingly, in the inset of Fig|2jb) are also drawn the 
pressure variations of two scaling quantities with respect 
to the Po = GPa state: s 
the boson peak and s^{P) 
bye frequency. As recently observed in [22| for a perma- 
nently densified Ge02 glass, it appears that the rescaled 
frequency of the boson peak s^(P) has a slope greater 
than the Debye one s^{P), which also points out that the 
CMT cannot fully account for the pressure dependence 
of low-vibrational dynamics in silica glass. 

Furthermore, it has been also shown in [T2j by studying 
the half-width Raman intensity main band, that amor- 
phous silica undergoes strong inter-tetrahedral structural 
rearrangements upon compression, which also shows that 
this glass is not homogeneous at nano-metric scale. We 
also measured the shift of the maximum of the inter- 
tetrahedra angle distributions with increasing pressure 
(not shown), and found values in agreement with exper- 
imental ones [l^ . 

In the inset of the FiglJl the ratio between the num- 
ber T+ of connected tetrahedral units that have the same 
rotation, and those t_ that have opposite ones, is plot- 
ted versus the hydrostatic pressure, after an elastic elon- 
gation has been applied to the material. This ratio 
T+ / T_ (P) is an indicator of the amount of inhomoge- 
neous atoms rearrangements at a given hydrostatic pres- 
sure state, namely when t_ > t+. In the inset of this 
figure, one can observes that the number of identical rota- 
tions T_|_ has a tendency to dominate the number of oppo- 



site inter-tetrahedral moves, when the hydrostatic pres- 
sure increases. In that way, the silica glass undergoes a 
decrease of inhomogeneous particle rearrangements with 
increasing pressure, i.e. the response of the glass is more 
and more homogeneous at all scales. This could explain 
the decrease of the boson peak intensity, as well as also 
introduces the possible link between structural effects re- 
lated to the disorder, and the vibrational properties of 
this glass. 

Finally, on the FiglUb), the vDOS of a Lennard- Jones 
glass under a null equivalent pressure density is also plot- 
ted, as well as the Debye prediction Eq.(|3]). The differ- 
ence that appears between the Debye's prediction and 
the computed vDOS is clearly drawn in the inset of the 
Figlljb), where the ratio between both quantities depicts 
a bump around ~ 1 rad.T~^ , to which is associated the 
boson peak anomaly. Also drawn on this FiglTJb) are the 
estimated contribution of longitudinal (L) and transverse 
(T) pulsations ojt/l{P) = 2TrCT/L{P)linaff{P) for the 
characteristic size ^„a/ / [P] of the elastic heterogeneous 
domains, already discussed in the previous Sec HI It then 
appears that the boson peak regime can be well approxi- 
mated by the lower bound pulsation u!t{P) of transverse 
waves associated with EH. 



B. 



Pressure dependence of the boson peak 
pulsation 



The last point discussed in the previous Sec lIII Al is 
more precisely approached on the Fig|4l^a) and (b). On 
these figures, for both glasses the boson peak pulsation 
wsp(P) measured on the Fig [1] is plotted in comparison 
to the pulsations ujt/l{P) for transverse and longitudi- 
nal elastic waves and for different pressures, i.e. using 
Cx/l{P)- Note that from previous studies, the charac- 
teristic sizes of EH for both glasses has been found to 
decrease and saturate at high pressure. On this figure, 
the power- law variation of the boson peak position ubp 
with P, as predicted by the soft potential model (SPM) 
j28l [29j , is also drawn. Our results seem to follow the 
expected P^/^ variation from SPM, thus assuming that 
bulk modulus arc pressure independent. For instance, 
in the case of silica, this situation is not fulfilled in the 
compression elastic regime. 

Then, on the same Fig|4j we see that the transverse 
pulsations trend, with increasing pressure, follows the bo- 
son peak one, while discrepancies appear for longitudinal 
pulsations. This emphasizes previous results j8l4lOl|: (i) 
the response of a glass to a long wavelength excitation is 
dominantly transverse and leads to the creation of noisy 
rotational atomic rearrangements, (ii) the characteristic 
size ^naff of these structures decreases and saturates at 
high pressure, (Hi) the boson peak is sited close to the 
resonant regime of transverse vibrational modes interact- 
ing with these structures or in other worlds of pulsations 
associated with their characteristic size. This last point 
is mainly phenomenological, but it bridges the scale be- 
tween the mechanical properties of glasses and their vi- 
brational spectrum. In the next section, we'll go deeper 
in this understanding by bridging the scale to the vibra- 
tional eigenvectors. 
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Fig 4. Variation of the measured boson peak pulsation 
ujbp{P) from Fig. [T] and for (a) 2D Lennard- Jones and (b) 
silica glasses, versus the hydrostatic pressure P. Boson peak 
pulsation is compared for each pressure to the estimation of 
the transverse ujt{P) and longitudinal ujl{P) pulsations asso- 
ciated with the elastic heterogeneities (EH) of characteristic 
size ^naff{P)- Note that from previous studies, the £,naff{P) 
decrease with increasing pressure. The P-dependence predic- 
tion of Gurevich et al. [28 . [2^ is also drawn, which quantita- 
tively scales with u)bp{P) and ujt[P)- 



IV. RESULTS FOR THE VIBRATIONAL 
EIGENVECTORS IN 2D 



«Rp(p) 



COrp(P) 




Fig 5. Participation ratio Prito, p) from Eq.lO, for two LJ2D 
systems under different densities. The corresponding pressure 
states are p = 0.9 cr"^ <^ P ~ Ocr^/e and p = 1.4 cr"^ P ~ 
130(7^/e. Only the first 2000 eigenvectors were used for this 
plot, and the boson peak pulsation uibp{p) from Fig[ljb) is 
also plotted. The change in the participation regime trend is 
marked by the boson peak pulsation, measured on the vDOS. 



seems to be more or less located in a region close to the 
lofFe-Regel crossover [23 . [sol WK . [38| , with wave excita- 
tions of predominantly transverse character. Recently, 
numerical studies (ill l30l | have been employed to confirm 
this last point, by computing the dynamic longitudinal 
and transverse structure factors (only the former is ac- 
cessible experimentally in the studied frequency range) 
from the knowledge of particle trajectories. In [ll[, the 
authors also showed using per particle vDOS, that trans- 
verse modes involved in the boson peak are associated 
with defective structures in the materials. 

In order to connect previous works related with EH 
in two-dimensional Lennard- Jones glasses (LJ2D) [sl-fioj 
and their characteristic size ^„a/ / (P) , the problematic is 
here approached by the use of the eigenvectors and their 
interaction with such specific zones. 

To this aim, we first analyse the participation ratio of 
eigenvectors along the vibrational spectrum. This one is 
obtained by computing the quantity: 



As discussed in Sec lIIBl the vibrational eigenvectors 
in the framework of the continuous elasticity theory are 
plane waves. In a disordered glassy material, one can ex- 
pect plane waves in the highest wavelength limit, while 
things may be more complicated with the decrease of 
the wavelength. Nevertheless, recent progress in com- 
puter simulations and experimental techniques have pro- 
vided new insights in the understanding of this prob- 
lem. For instance, in silica glasses j33l - l35l |. it has been 
shown that acoustic phonon waves undergo a crossover 
between a propagation regime and another in which they 
are strongly diffracted. This crossover corresponds to 
the loffe-Regcl limit where the phonons mean free path 
is of same order than their wavelength [s^, [s^ , while the 
boson peak has been found to be closely related to this 
crossover frequency. It then appears that the boson peak 



-l E^^ I e^(a;,p) P 
Ef |e.(a.,p) 1^ 



Pr{,UJ,p)^N 



(5) 



where e{uj,p) are the eigenvectors for a pulsation mode 
Lo under an overall hydrostatic pressure P or density p. 
When a mode is delocalized and then most of atoms vi- 
brate to the corresponding pulsation a;, Pr{uj,p) ~ 1, 
however, if the mode is localized, then few atoms par- 
ticipate to the mode and Pr{uJ,p) — > N~^. The par- 
ticipation ratio is plotted on the FiglS] for two density 
values corresponding to a free p = 0.9 cr~^ P ^ Oa^ /e 
state and a strongly compressed one with p = lAa^^ <^ 
P ^ 130a^/e. On this figure, it appears that Pr{uj,p) 
decreases from a full atomic vibrational participative 
state to a lower one that is stabilized around 30% over 
the studied pulsation range. Let us remind that for 
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Fig 6. Noise in eigenvectors from Eq.((7]) versus pulsation uj 
for different glass densities in 2D. For each density, the mea- 
sured boson peak pulsation ujbp{p) is marked with a vertical 
arrow, while the vertical dashed line is the value of the pulsa- 
tion liJt[p) = 2-kCt[p) / inaf f{p) that would have a transverse 
wave of wavelength of the same size than £,naff{p)- The noise 
becomes predominant for pulsations close to ujt{p) and sat- 
urates when the boson peak ubp{p) is reached. Inset: pro- 
jected amplitudes over plane waves from Eq.®, and for a 
null equivalent hydrostatic pressure. The distance to the the- 
oretical plane wave is given by {pth — p)'- p = 3 is the first 
mode after the two Goldstone's ones, p = 43 is somewhere 
before the boson peak, and p = 600 is far after the boson 
peak. Spreading of the projected amplitude over neighbour- 
ing modes is observed, which increases with the considered 
mode index p. For the highest mode p = 600, the plane wave 
decomposition does not hold. 

the P ^ Ocr^/e sample, the Debye pulsation is around 
(jjjj{P) Pi Ih.Zrad.T^^ . Then, one may expect for higher 
pulsations Pr(w, p) to tend to A^"^. A crossover appear in 
Pr at a pulsation close to the boson peak one (measured 
using the vDOS) which means that around the boson 
peak, a clearly non-null amount of atoms participate to 
the vibrational states. 

A closer look to the eigenvectors can be obtained 
by defining a pseudo-spectral density for an eigenstate 
e(a;p, r, p)) for a mode p of pulsation uip, over a theoret- 
ical plane wave of wave- vector kg, given for a mode q of 
pulsation w^: 

Ap{q) = (e(tjp,r,p) | etheo{i^q,r, p)) (6) 

With the above definition, one can build the noisy eigen- 
vector (5e(ajp, r, p): 

5e{ujp, r, p) = e(wp, r, p) - ^ Ap{q)etheo{t^q, r, p) (7) 

where T>p is the degeneracy set for the mode p of pulsa- 
tion u)p. On inset of the Fig|6l the pseudo-spectral den- 
sity Eq.® is plotted for the P ^ Ocr^/e sample, versus 
the "distance" [ptu — p) to the theoretical mode q — pth, 
and for three modes: p = 3 is the first after the two first 
Goldstone's modes, p = 43 is an arbitrary mode before 



Fig 7. Relative noise in eigenvectors as a function of the 
reduce wavelength A(a;,p) (see text for details), and for 
two samples under different density. The Rayleigh regime 
(~ A"'^), where plane waves are diffracted by the weak dis- 
order, holds for wavelengths greater than EH size. An inter- 
mediate exponential regime forms a crossover for wavelengths 
between EH size and boson peak wavelength, while strong 
diffraction occur for lower wavelengths. 



the boson peak, and p = 600 is far away after the boson 
peak mode. In the inset of this Figl6l one can note that 
the lowest mode amplitude is well stuck over the theo- 
retical plane wave, while increasing the mode index leads 
to a spreading of the projected amplitudes Ap{pth —p)- 
This means that the plane wave picture used to describe 
eigenvectors is more and more lost when the mode index 
increases, or at least that as soon as the mode index in- 
creases, more and more neighbouring waves pollute the 
ideal theoretical plane wave which should represent the 
p mode. This can be more clearly seen by computing the 
noise in eigenvectors from Eq.©. Results are shown on 
the main panel of the Fig|6]for three different hydrostatic 
pressures. The trend of 6e(LUp, r, p) is compared for each 
density p to the boson peak pulsation ub p (p) , and to the 
pulsation wt(p) = '^'^C^xip) / £,naf f (p) of resonant trans- 
verse waves of wavelength equal to the EH characteristic 
size ^.naffip)- It then appears that the noise becomes 
predominant for pulsations around uit{p) and that start- 
ing from the boson peak regime, the plane wave picture 
no more holds [H |M HI iH^ . 

An interesting point also results from the quantity 
ujt{p)- If for all considered densities, this one lies just 
before the noisy saturated regime, one can also wonder 
how the EH affect the plane wave framework and thus 
are involved in the propagation of elastic waves. On the 
FigEl the problematic of the interaction process of plane 
waves with EH is then approached by studying the noise 
relative ratio || Se{uj,r,p) \\ / \\ e(w,r, p) ||. On this fig- 
ure, a reduced wavelength A(u!,p) = Xtheoi'-^, p)/^l,t{p) 
is used, with ^l,t{p) the scaling lengths from the 2D 
scaling mode analysis @ for longitudinal (L) and trans- 
verse (T) waves, and Xtheo{'-^iP) the theoretical wave- 
length for a plane wave at a given pulsation uj in the con- 
tinuum medium under a density p. Note that we have 
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observed that the EH characteristic size £,naff{p) — £,t{p) 
while £,l{p) ^ 2^t(p), and that LO*j^{p) — oJxip) with 
^L.rip) " '^'^Cl,t{p)/^l,t{p)- Hence, with such a choice 
of units, one can remove the longitudinal and trans- 
verse wave types under an unique plot, while keeping 
in mind that one can recover true units by using the 
relation between £,naff{p), and ~ 2^t{p)- On 

the same figure, the reduce wavelength for the boson 
peak A{ujbp,p) is also drawn, with the true wavelength 
\{ujbp,p) = '2t^Ct,l/£,t,l{p)- Again, using the reduced 
wavelength A(aj,p), the boson peak wavelength thereby 
estimated is the same for (L, T) waves. On this Figl?! 
it appears three regimes for the relative noise amplitude: 
(i) for wavelengths greater than the EH characteristic 
size, a Rayleigh type diffusion process is observed, i.e. 
planes waves are scattered by a weak disorder, (ii) for 
wavelengths between the boson peak one and ones of the 
same size than EH, an intermediate exponential regime 
seems to occur, (in) for wavelengths lower than the boson 
peak one, the noise is of same order than the eigenvector 
and plane waves are strongly scattered by the disorder. 
The same scenario holds for an higher pressure sample, 
as shown in the inset of the Fig 17] 

A useful way to characterize more precisely the inter- 
action process between plane waves and disorder, is to 
define a scalar correlation function of the noisy eigenvec- 
tors given in Eq.©. This one is then written as: 




A(co,p) 



10"' lo" lO' 

A(CO,p) = \heo^®'P^ ' ^T,L^P) 



Fig 8. Dominant rotational structure size in noisy eigenvec- 
tors 5e{ijj,r,p) from the first zero of the scalar correlation 
function Eq.®, versus the reduced wavelength K{uj,p) for 
T and L waves (remind that ~ 2^t). (C^(p))t^ is equiv- 
alent to the threshold mean free path length for propagat- 
ing phonons in disordered media. Three regimes occur with 
crossovers between each marked by, first EH characteristic size 
at A.{u),p) — 1, second the boson peak position A(cj_Bp,p). 
The boson peak is sited at the loffe-Regel [s^ limit. In- 
set: same plot for system under an equivalent pressure of 
P ~ 130cr/e^ 



CUr,p) = {Se{uj,r,p)\Se{uj,0,p)) 



(8) 



This correlation function depicts an anti-correlation 
zone at a distance r = (C^(p))cj where the average is 
taken over the degeneracy set Vp for modes p of same 
pulsation ujp = lu. This crossover from a positive cor- 
relation to a negative one, which thus define the first 
zero {C^{p)) uj, is a track of the emergence of a dominant 
rotational structure in the noisy eigenvectors de{uj,r, p), 
and {C^{p))u: can then be considered as a characteristic 
or threshold mean free path length for phonons propaga- 
tion. 

On the Figig the length (CHp))'^ is plotted 0jS 0, func- 
tion of the reduced wavelength A(aj, p) for both longitu- 
dinal (L) and transverse (T) waves, and for two differ- 
ent sample density. On the same plot, the boson peak 
position A{iOBP,p) is drawn. Again, three regimes are 
observed: (i) for phonon wavelengths greater than the 
characteristic length of EH at A{u!,p) = 1, the thresh- 
old mean free path length is independent of the mode, 
meaning that plane waves propagation is not strongly 
altered by the disorder; (ii) A(lu,p) = 1 marks an inde- 
pendent/dependent crossover until A — >■ A{ujbp,p), and 
during which a linear dependence of the threshold mean 
free path length {C^{p))u: with the wavelen gth appears: 
this is equivalent to the loffe-Regel regime 
(Hi) for wavelengths lower than the boson peak one, the 
strong diffraction regime of plane waves by the disor- 
der is reached, with a quadratic dependence of the mean 
free path length with phonon wavelength. The boson 
peak then marks a crossover at the loffe-Regel limit 
tni [13, in, m [13. Note that the same scenario also 
holds for a more compressed sample as shown in the in- 



set of the FiglHl 

Finally, on the FigO results are summarized for dif- 
ferent quantities encountered during this work. On this 
figure, data are plotted as the inverse of a specific length, 
namely by considering the ratio {lux /Ct){p) oc l^^{p). 
As previously discussed in [l^, [lH H^l , and in agreement 
with the present work, it appears that transverse modes 
are more affected by the disorder, as well as more involved 
in the boson peak anomaly than longitudinal ones. Con- 
sequently, this motivates the choice of transverse speed 
of sound Ct{p) to perform the above ratio. The quan- 
tities X are related to: (a) EH characteristic size and 
X := {T,naff), (b) boson peak and X := (BP), (c) 
the threshold mean free path (C^(p))w in the mode inde- 
pendent regime, i.e. for wavelengths larger than ^na/ / (p) , 
and X := (T, vec). On the Fig[51 we also plotted the in- 
verse boson peak wavelength computed from the trans- 
verse power spectrum: 



N 



2^kAe(a;,rj 



p)e 



ik.r , 



(9) 



where k = k/ || k ||. Then the transverse boson peak 
wavelength is computed by MAXk{/T(wBp, k, p)}, giv- 
ing II ksp 11= 2tt/Xt,bp[p)- Hence, on the Fig[Sl the 
pressure dependence of all these quantities are drawn. 
Interestingly, the inverse boson peak length Xt,bp{p) 
extracted from Eq.([5|) seems to not be altered by the 
pressure variation, and is very close to the inverse boson 
peak length computed using transverse speeds of elastic 
waves. The comparison between both independent mca- 
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Fig 9. Plot of the inverse lengths from (ujx /Ct){p) ratio. De- 
pending X observable, lengths represent: EH characteristic 
size with uJT,naff, boson peak pulsation ujbp{p), lowest mode 
threshold mean free path for elastic waves with uJT,vec- Main 
panel: for the LJ glass, comparison is also shown with the 
boson peak wavelength, independently computed at a; = ujbp 
from the transverse eigenvectors power spectrum given in 
Eq.(|9l). The weak variation with varying pressure points out 
the fact that the boson peak is a signature of a characteristic 
length which emerges from an ordering/frustration competi- 
tion in glasses. Inset: same plot for the silica glass. 



surements seems to point out the fact that the pulsation 
shift of the boson peak is directly related to the increase 
of speed of elastic waves under compression. We then also 
point out that the same trend seems to occur for inverse 
lengths computed from EH characteristic size and lowest 
mode threshold mean free path for elastic waves. In the 
inset of the FiglHl the same plot is drawn for the silica 
glass, and the same trend can also be observed. This 
type of plot should then encourage to view the boson 
peak anomaly as a process that results from the exis- 
tence of a characteristic length, itself inherent from the 
competition between ordering and frustration in glassy 
materials. 



V. CONCLUDING REMARKS 

During this work, we have extended to the pressure 
dependent regime, previous studies on a prototypical 
two-dimensional soft glass former. Taking advantage 
of known results concerning its mechanical response to 
a macroscopic excitation, and the relation it exists be- 
tween the involved specific non-affine atomic displace- 
ments (which are connected with the existence of EH) 
and its vibrational spectrum, we have developed an ad- 
ditional point of view on the boson peak scenario. This 
one has been achieved by studying the interplay between 
the vibrational eigenvectors and the well characterized 
EH, for this kind of system under continuous compres- 
sion. 

We have demonstrated that eigenvectors with wave- 
lengths larger than the EH characteristic size can be well 



approximated by plane waves, as predicted by the con- 
tinuum theory of elasticity. Conversely, when the wave- 
length becomes of same order than this characteristic 
size, a crossover occurs and the plane wave framework no 
more holds. Hence, elastic waves are strongly scattered 
by the EH, while this crossover is marked by the boson 
peak pulsation, that also coincides with the lofFe-Regel 
limit. 

We have also shown that the boson peak pulsation 
shifts with the pressure on a same way than the pulsa- 
tion of transverse elastic waves of wavelength of same size 
than EH, for both soft LJ and strong silica glasses. This 
results confirms the fact that, whatever the pressure, the 
boson peak is always sited close to the non-affine regime, 
where the response of the system to a large wavelength 
excitation is self-organized in a collective manner and in 
such a way that it emerges a characteristic length as the 
result of the competition between order and frustration 
in glasses. Then, we also pointed out that the pulsation 
shift of the boson peak is directly related to the increase 
of the speed of sound for transverse elastic waves. It fol- 
lows that the boson peak can be indeed viewed as the 
result of a length, which might be related to this emerg- 
ing length. 

A part of this work was also devoted on the preliminary 
study of the pressure dependence of the eigcn-frequcncies 
of a silica glass, a strong glass former. The density of 
state of such a glass was computed for various pressures, 
while the mechanical response of this glass to a macro- 
scopic elongation was also performed in a parallel inde- 
pendent study, in order to extract the characteristic size 
of EH. We have then measured a decrease and a satura- 
tion of this characteristic length with increasing pressure. 
By the same way, it has been also possible to extract the 
inter-tetrahedral reciprocal rotations in response to the 
macroscopic deformation. The number of identical ro- 
tations has been shown to dominate the opposite ones 
when the pressure increases, which can be explained by 
a decrease of inhomogeneous moves, on the same manner 
than EH do. The decrease of the boson peak intensity 
and its shift are usually explained by a suppression of 
modes responsible for the boson peak. Another expla- 
nation could be based on the variation of inhomogeneous 
moves and EH when the pressure increases, as in the case 
of the soft LJ glass former. 

Hence, compiling these results with the ones obtained 
from the vibrational spectrum, it has been shown that 
the boson peak frequency follows the same pressure de- 
pendence than that of transverse elastic waves with fre- 
quency associated with the EH characteristic size mea- 
sured in that glass. This corroborates results from the 
simple two-dimensional soft LJ glass former, and should 
encourage to extend the methodology developed for the 
eigenvectors study, to the case of the silica glass. The 
task is much more harder computationally, as one has 
to consider systems with a large amount of particles in 
order to get rid of finite size effects . It should be also 
interesting to have a more precise look at the involved 
structural changes when the pressure increases, and how 
they could be related with elastic heterogeneities. 
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